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Spherically symmetric solutions in a FRW background 
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We impose perfect fluid concept along with slow expansion approximation to derive new solutions 
which, considering non-static spherically symmetric metrics, can be treated as Black Holes. We will 
refer to these solutions as Quasi Black Holes. Mathematical and physical features such as Killing 
vectors, singularities, and mass have been studied. Their horizons and thermodynamic properties 
have also been investigated. In addition, relationship with other related works (including mcVittie’s) 
are described. 


I. INTRODUCTION 


The Universe expansion can be modeled by the so called FRW metric 

dr 2 

ds 2 = —dt 2 + a(t) 2 [- -—+ r 2 d0 2 + r 2 sin{9) 2 d(j) 2 ], (1) 

(1 — kr~) 

where k = 0,+1,—1 are curvature scalars which represent the flat, closed and open universes, respectively. The 
WMAP data confirms a flat (k = 0) universe JjJ. a(t) is the scale factor and for a background which is filled by a 
perfect fluid with equation of state p = ojp , there are three classes of expanding solutions. These three solutions are 


a(t) = oof 3( ‘ J+1) 


for oj ^ 0 when — 1 < oj and , 


a{t) = cioe Ht 


for oj = —1 (dark energy), and for the Phantom regime (oj < —1) is 

a(t ) = oq (to - t) 3( “ +1) , 


( 2 ) 

( 3 ) 

( 4 ) 


where to is the big rip singularity time and will be available, if the universe is in the phantom regime. 

In Eq. ([3]), H(= ^y) is the Hubble parameter and the current estimates are H = 73t_\kms~ l Mpc~ x [l.|. 

Note that, at the end of the Phantom regime, everything will decompose into its fundamental constituents Q. In 
addition, this spacetime can be classified as a subgroup of the Godel-type spacetime with a = m = 0 and k' = 1 [3]. 

A signal which was emitted at the time to by a co-moving source and absorbed by a co-moving observer at a later 
time t is affected by a redsliift (z) as 


1 + z = 


a(t) 

a(t 0 )' 


( 5 ) 


The apparent horizon as a marginally trapped surface, is defined as [!] 


= o, 


which for the physical radius of £ = a(t)r , the solution will be: 




1 



( 6 ) 

( 7 ) 
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The surface gravity of the apparent horizon can be evaluated by: 

n = =d a (\T-hh ab d b i ;). 


( 8 ) 


Where the two dimensional induced metric is h a b = diag(— 1, It was shown that the first law of thermo¬ 
dynamics is satisfied on the apparent horizon The special case of w = — 1 is called the dark energy, and by a 

suitable change of variables one can rewrite this case in the static form jj|: 

ds 2 = -(1 - H 2 r 2 )dt 2 + ——+ r 2 dfl 2 . (9) 

(1 — H-r-) 

This metric belongs to a more general class of spherically symmetric, static metrics. For these class of spherically 
symmetric static metrics, the line element can be written in the form of: 

dr 2 

ds 2 = -f(r)dt 2 + ttt + r 2 dfl 2 , (10) 

f{r) 

where the general form of f(r) is: 

/(r) = l-2^|-ffV. (11) 


In the above expression, m and Q represent mass and charge, respectively. For this metric, one can evaluate redshift: 


1 + 2 


1 _ 2™ + 91 _ H 2 r 2 i 

/ _ r _ 7 __\ 7 

k T-rn O ' 


1 - 2 - 


+ ^ — H 2 r 2 ' 


( 12 ) 


Where, ro and r are radial coordinates at the emission and the absorption points. For the horizons, the radius and 
the surface gravity can be found using equations 

gtt = 

k = 

where (') denotes derivative with respect to the coordinate r [9|. From the thermodynamic laws of Black Holes (BHs) 
we know 


f[r) =0 —> Th (13) 

f'(r) 

2 Wh ’ 



(14) 


which T is the temperature on the horizon Q • Validity of the first law of the thermodynamics on the static horizons 
for the static spherically symmetric spacetime has been shown 0,0. 

The BHs with the FRW dynamic background has motivated many investigations. The first approach, which is 
named Swiss Cheese, includes efforts in order to find the effects of the expansion of the Universe on the gravitational 
field of the stars fl2j . introduced originally by Einstein and Straus (1945) [l3j]. In these models, authors tried to 
join the Schwarzschild metric to the FRW metric by satisfying the junction conditions on the boundary, which is 
an expanding timelike hypersurface. The inner spacetime is described by the Schwarzschild metric, while the FRW 
metric explains the outer spacetime. These models don’t contain dynamical BHs, Because the inner spacetime is in 
the Schwarzschild coordinate, hence, is static 0- In addition, the Swiss Cheese models can be classified as a subclass 
of inhomogeneous Lemabitre-Tolman-Bondi models [THI . fldj . 

Looking for dynamical BHs, some authors used the conformal transformation of the Schwarzschild BH, where 
the conformal factor is the scale factor of the famous FRW model. Originally, Thakurta (1981) have used this 
technique and obtained a dynamical version of the Schwarzschild BH 0- Since the Thakurta spacetime is a conformal 
transformation of the Schwarzschild metric, it is now accepted that its redshift radii points to the co-moving radii of 
the event horizon of BH 0,0,0|. By considering asymptotic behavior of the gravitational lagrangian (Ricci scalar), 
one can classify the Thakurta BH and its extension to the charged BH into the same class of solutions 00. The 
Thakurta spacetime sustains an inward flow, which leads to an increase in the mass of BH fl8l - |20i ] . This ingoing flow 
comes from the back-reaction effect and can be neglected in a low density background (20|. In fact, for the low density 
background, the mass will be decreased in the Phantom regime fut . Also, the radius of event horizon increases with 
the scale factor when its temperature decreases by the inverse of scale factor 00. 
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Using the Eddington-Finkelstein form of the Schwarzschild metric and the conformal transformation, Sultana and 
Dyer (2005) have constructed their metric and studied its properties (22|. In addition, unlike the Thakurta spacetime, 
the curvature scalars do not diverge at the redshift singularity radii (event horizon) of the Sultana and Dyer spacetimes. 
Since the Sultana and Dyer spacetimes is conformal transformation of the Schwarzschild metric, it is now accepted 
that the Sultana and Dyer spacetimes include dynamic BHs [3. Various examples can be found in [n| l23l425j | . 
Am ong these conformal BHs, only the solutions by M c Clure et al. and Thakurta can satisfy the energy conditions 


16J19L Static charged BHs which are confined into the FRW spacetime and the dynamic, charged BHs were studied 
2ff - 33] . The Brane solutions can be found in 


In another approach, mcVittie found new solutions including contracting BHs in the coordinates co-moving with 
the universe’s expansion f37|. Its generalization to the arbitrary dimensions and to the charged BHs can be found in 
[38 . 39). In these solutions, it is easy to check that the curvature scalars diverge at the redshift singularities. In this 
approach, authors have used the isotropic form of the FRW metric along as the perfect fluid concept and could find 
their solutions which can contain BHs J40|. The mass and the charge of their BHs seem to be decreased with the scale 
factor. Also, it seems that the redshift singularities does not point to a d yna mic event horizon frO - tid] ]. Unlike the 
Swiss Cheese models, the energy conditions are violated by these solutions [lfj. These solutions can be considered as 
Models for cosmological inhomogeneities fl5j . 

This paper is organized as follows: in the next section, we consider the conformal transformation of a non-static 
spherically symmetric metric, where conformal factor has only time dependency. In addition, we derive the general 
possible form of metric by using perfect fluid concept. In section 3, slow time varying approximation is used in order to 
find the physical meaning of the parameters of metric. In continue, the mcVittie like solution and its thermodynamic 
properties are addressed. In section 4, we generalized our debates to the charged spacetime, when the effects of the 
dark energy are considerable. In section 5, we summarize and conclude the results. 


II. METRIC, GENERAL PROPERTIES AND BASIC ASSUMPTIONS 


Let us begin with this metric: 

dr 2 

ds 2 = a(r) 2 [-/(t, r)dr 2 + T - . -r + r 2 dd 2 + r 2 sin(6) 2 dc/) 2 }. (15) 

(1 - kr 2 )f{r,r) 

Where a(r) is the arbitrary function of time coordinate r. This metric has three Killing vectors 

<9^, sin cj) dg + cot 9 cos ft and cos <j> dg — cot 6 sincj) d^. (16) 

Now, if we define new time coordinate as 


T — } t — 


a(r)dT, 


(17) 


we will get 

„ dr 2 

ds 2 = -f(t, r)dt + a(t) 2 [ _ kr 2 )f(t r) + + r sin ( e ) ^ ^ ( 18 ) 

which possesses symmetries like as Eq. ©. From now, it is assumed that a(t) is the cosmic scale factor similar to 
the FRW’s. For f(t,r ) = 1, Eq. (fTBl) is reduced to the FRW metric ([]]). Also, conformal BHs can be achieved by 
choosing f(t,r) = f(r) where, the general form of /(?') is fu| : 


/(r) = 1 - 


2 m Q 2 


r — 


Ar 2 

— 


Therefore, conformal BHs can be classified as a special subclass of metric (fT51) . n a 
r = const and yields 


n a n a 


= 9 


(1 - kr 2 )f(t,r) 
a(t) 2 


(19) 

5 r a is normal to the hypersurface 


( 20 ) 


which is timelike when (1 — kr 2 )f(t, r) < 0, null for (1 — fc? ,2 )/(t, r) = 0 and spacelike if we have (1 — kr 2 )f(t, r) > 0. 
For an emitted signal at the coordinates to and r o, when it is absorbed at coordinates f and r simple calculations lead 
to 


l + z 


X_ 

Ao 


a(t) , f(t,r) i 
a(toU f(t Q ,r 0 ) 


( 21 ) 
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as induced redshift due to the universe expansion and factor f(t,r). Redshift will diverge when /(foUo) goes to zero 
or 1 + z —> oo. This divergence as the signal of sing ularity is independent of the curvature scalar (fc), unlike the 
Mcvittie’s solution and its various generalizations 38, ;39], which shows that our solutions are compatible with the 
FRW background. As a desired expectation, it is obvious that the FRW result is covered when /(to, r o ) = f{t, r) = 1. 
The only non-diagonal term of the Einstein tensor is 


G tr = — 


1 — kr 2 
f(t,r)a{t) 3 r 


( a{t)f(t,r) - f'(t,r)a(t)r), 


( 22 ) 


which (') and (') are derivatives with respect to time and radius, respectively. Using one gets 


G tr = - h:1]n{t l (a(t)f(a(t),r) - 


f(a(t),r)a(t) 3 r 


(23) 


where f(a(t),r) = In order to get perfect fluid solutions, we impose condition G tr = 0 and reach to 


= f(a(t)r) = ^ b n (a(t)r) n . 


(24) 


Although Eq. (El) includes numerous terms, but the slow expansion approximation helps us to attribute physical 
meaning to the certain coefficients b n . Since Gt r = 0, we should stress that here that there is no redial flow and 
thus, the backreaction effect is zero [l)| [20[, which means that there is no energy accretion in these solutions [45| . 
Finally and briefly, we see that the perfect fluid concept is in line with the no energy accretion condition. The only 
answer which is independent of the rate of expansion can be obtained by condition b n = S n o which is yielding the 
FRW solution. 


III. MCVITTIE LIKE SOLUTION IN THE FRW BACKGROUND 


The mcVittie’s solution in the flat FRW background can be written as [16) 


ds 2 = -( 


1 - 

2 a(t)r 

1 + 

M 


M 


) 2 dt 2 + a(t) 2 ( 1 + ^^ndf 2 +f 2 dfl 2 ]. 


This metric possess symmetries same as metric m- f is isotropic radius defined by: 

AT 


r = f:(1+ 2f) ■ 


(25) 


(26) 


There is a redshift singularity at radii fh = ^tj which yields the radius = ^^r(l + a{t)) 2 [51]- In addition, fh is 
a spacelike hypersurface, and can not point to an event horizon j4f|. 

Consider f(a(t)r) = 1 — This assumption satisfies condition (El) and leads to 


d » 2 =-(* - wt ^ ++ ^ (27) 

For b -1 0, this metric will converge to the FRW metric when r —> oo. The Scliwarzschild metric is obtainable by 

putting a(t) = 1, b-i = M and k = 0. Metric suffers from three singularities at a(t) = 0 (big bang), r = 0 and 


f{a(t)r) = 0 => a(t)r h = 26_i. (28) 

Third singularity exists if 6_i > 0. In this manner, Eq. (1211) will diverge at ro = r^. In addition and in contrast to the 
Gao’s solutions, the radii of the redshift singularity (r^) in our solutions is independent of the background curvature 
(k). while for the flat case our radius is compatible with the previous works [3, H3, S]. Also, metric changes its 
sign at r = rh just the same as the schwarzschild spacetime. In addition, curvature scalars diverge at this radius as 
well as the mcVittie spacetime. Accordingly, this singularity point to a naked singularity which can be considered 
as alternatives for BHs (3, E3]- In continue, we will point to the some physical and mathematical properties of this 
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singularity which has the same behaviors as event horizon if one considers slow expansion approximation. The surface 
area integration at this radius leads to 


A = J yfadOdcj) = 47rr^a(f) 2 = 167 t(&_i) 2 


(29) 


The main questions that arise here are: what is the nature of 6 _i? and can we better clarify the meaning of r^? For 

these purposes, we consider the slow expansion approximation (a(t) « c), define new coordinate 77 = cr and get 


ds 2 


— (1 — ———)dt 2 
V 


+ 


d?y 2 


(i-fcV)(i-^) 


r] 2 dd 2 + 77 2 sin(9) 2 d(j) 2 


(30) 


where k! = ^. In these new coordinates, (f, 77, 9, </>), and from Eq. (l20ll it is apparent that for 6_ 1 > 0, hypersurface 
with equation 77 = r]h = 26_i is a null hypersurface. When our approximation is broken, then r]h may not be actually 
a null hypersurface, despite its resemblance to that. We call this null hypersurface a quasi event horizon which is 
signalling us an object like a BH and we refer to that as a quasi BH. From now, we assume b- 1 > 0, the reason of this 
option will be more clear later, when we debate mass. Therefore by the slow expansion approximation, (= 2b ~ 1 ) 
plays the role of the co-moving radius of event horizon and it is decreased with time. In order to find an answer to 
the first question about the physical meaning of &_i, we use Komar mass: 


M = f n a crpVat; 1 / 1 dA, 
J s 


(31) 


where is the timelike Killing vector of spacetime. Since the Komar mass is only definable for the stationary and 
asymptotically flat spacetimes 0 ], one should consider the flat case (k = 0 ) and then by bearing the spirit of the 
stationary limit in mind (the slow expansion approximation) tries to evaluate Eq. (1311) . 

Consider n a = — |and 073 = , ^ as the unit timelike and unit spacelike four-vectors, respectively. 

Now using Eq. m and bearing the spirit of the slow expansion approximation in mind, one gets 


M= 4~ n °‘ a ( )T <*t d A = 6 -i> 
J s 


(32) 


which is compatible with the no energy accretion condition ( Gt r = 0). In addition, we will find the same result as 

1 


-irwgfy- Since 

V 1 v 


Eq. (13^1) . if we considered the flat case (k = 0) of metric (l30l) and use n a = 1 — —^<5^ and 073 = 

the integrand is independent of the scale factor (a(t)), the slow expansion approximation does not change the result of 
integral. But, the accessibility of the slow expansion approximation is necessary if one wants to evaluate the Komar 
mass for dynamical spacetimes J48|. Indeed, this situation is the same as what we have in the quasi-equilibrium 
thermodynamical systems, where the accessibility of the quasi-equilibrium condition lets us use the equilibrium for¬ 
mulation for the vast thermodynamical systems [49[. It is obvious that for avoiding negative mass, we should have 
b -1 > 0. Relation to the Komar mass of the mcVittie’s solution can be written as jla. 1 391 


M, 


mcVittie — 


M 

a {t)' 


(33) 


In addition, some studies show that the Komar mass is just a metric parameter in the mcVittie spacetime mmm- 
Indeed, Hawking-Hayward quasi-local mass satisfies M = 0, which is compatible with G r t = 0 and indicates that 
there is no redial flow and thus the backreaction effect, in the mcVittie’s solution (l9l - !2ll |45| . In order to clarify 
the mass notion in the mcVittie spacetime, we consider the slow expansion approximation of the mcVittie spacetime 
which yields 


ds 2 


1 _ M 

-( 2r > ) 
'1 1 M _ r 


2 dt 2 + (1 + ^) 4 [dT 7 2 + rj 2 dn 2 ]. 
2?7 


(34) 


This metric is signalling us that the M may play the role of the mass in the mcVittie spacetime. In addition, by 
defining new radii R as 

R{t,r) = a(t)f( 1 + ^) 2 , 


( 35 ) 
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one can rewrite the mcVittie spacetime in the form of 

ds 2 = _(i _ ™ _ H 2 R 2 )dt 2 - 2HR dtdR + + ^ 2 dn 2 , 

* 1 


( 36 ) 


where H = - [HO]. This form of the mcVittie spacetime indicates these facts that the Komar mass is a metric 
parameter and M is the physical mass in this spacetime [50j |. Finally, we see that the results of the slow expansion 
approximation (Eq. (13411 ) and E g. (13611 are in line with the result of the study of the Hawking-Hayward quasi-local 
mass in the mcVittie spacetime [44], |42j, (45|, [5C| ■ For the flat case (k = 0) of our spacetime (Eq. (1271) 1. by considering 
Eq. (1331) and following the slow expansion approximation, we reach at 

ds 2 w —(1 — ‘^-)dt 2 H- ^ 2 M 7 7 2 ^ 2 + 7 f sin(9) 2 deft 2 . (37) 

^ (1 r? ) 

Also, if we define new radius R as 


we obtain 


R M . 2 

r = V 1+ 2«)' 


(38) 


ds 2 


,(i zMl _ m 2 ( i + ^) 6 v ,,2 _ 2i?g(i + #) 5 

(1 + ^) 2 ( I -#) 2 ( 1 - H ) 

(1 + ^) 4 [dR 2 + R 2 dn 2 }. 


(39) 


Both of the equations ( 1371 ) and ( 1551 ) as well as the no energy accretion condition suggest that, unlike the mcVittie’s 
spacetime, the Komar mass may play the role of the mass in our solution. From Eq. (1391) it is apparent that R = ^- 
points to the spacelike hypersurface where, in the metric (l36l) . R = 2 M points to the null hypersurface. In the 
next subsection and when we debate thermodynamics, we will derive the same result for the mass notion in our 
spacetime.Only in the a(t) = 1 limit (the Schwarzscliild limit), Eqs. (l39l) and (l25l) will be compatible which shows 
that our spacetime is different with the mcVittie’s. Let us note that the obtained metric (Eq. ( 15 ( 11 ) ) is consistent with 
Eq. (l36l) . provided we take M = 0 (the FRW limit). 


Horizons, energy and thermodynamics 

There is an apparent horizon in accordance with the FRW background which can be evaluated from Eq. ©: 

( 1 - fcr ap)( 1 - ) 2 ~ r 2 ap a{t) 2 = 0. (40) 

Oj\i)T ap 

This equation covers the FRW results in the limit of M —> 0 ( see Eq. J7J)). In addition, one can get the Schwarzscliild 
radius by considering a{t) = 0, which supports our previous definition for b-\. Calculations for the flat case yield four 
solutions. The only solution which is in full agreement with the limiting situation of the FRW metric (in the limit of 
zero M) is 


Tap = 


l + Vl^8HM 

2d 


Therefore, the physical radius of apparent horizon (£ op = a(t)r ap ) is 


(41) 




ap — 


1 + Vl - 8 HM 
2 H 


(42) 


which is similar to the conformal BHs Hi- It is obvious that in the limit of M —> 0, the radius for the apparent 
horizon of the flat FRW is recovered. For the surface gravity of apparent horizon, one can use Eq. d5J) and gets: 


(1 


kfrw__ 

2M \ 2 
a(t)r a p ' 


M 1 


2 M > 
a{t)r ap * 


(< d(t ) + 2 



(43) 


K = 


(i 
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where h ab = diag{— \ M , r ap is the apparent horizon co-moving radius m and kfrw is the surface 

1 TFTyr 1 > 

gravity of the flat FRW manifold 


a(t ) 2 + a(t)a(t) 

KFRW= - 2a(t)a(t) 


(44) 


The schwarzschild limit (k = -jjg) is obtainable by inserting a(t) = 1 in Eq. (H51) . In the limiting case M —► 0, 
Eq. (H51) is reduced to the surface gravity of the flat FRW spacetime, as a desired result. The Misner-Sharp mass 
inside radius £ for this spherically symmetric spacetime is defined as [521 ]: 


M M s = \{l-h ab d a ^d b 0- (45) 

Because this definition does not yield true results in some theories such as the Brans-Dicke and scalar-tensor gravities, 
we are pointing to the Gong-Wang definition of mass [53]: 


M GW = ^(i + h ab d a td b t). 


(46) 


It is apparent that, for the apparent horizon, Eqs. (HU) and (1461) yield the same result as Mow = Mms = In 
the limit of M —» 0 , the FRW’s results are recovered and we reach to M G w = Mms = pV as a desired result Ca¬ 
using Eqs. m and m and taking the slow expansion approximation into account, we reach to M G w = Mms — M 
as the mass of quasi BH. Also, this result supports our previous guess about the Komar mass as the physical mass 
in our solution, and is in line with the result of Eqs. (l37l) and (l39l) . For the Mcvittie metric, Eqs. (l45j) and (l46l) yield 
M G w = Mms — x as the confined mass to radius £/,, = a{t)fh = 4f ■ Also, Eqs. (15^1) . (H5l) and (H7>1) leads to the same 
result in the Schwarzschild’s limit (M = M GW = M M g = M). For the flat background, using metric (13U1) . Eq. (fl3l) 
and inserting results into Eq. one gets 


T ~ 


8nM 


(47) 


for the temperature on the surface of quasi horizon. The same calculations yield similar results, as the temperature 
on the horizon of the Mcvittie’s solution. For the conformal Schwarzschild BH, the same analysis leads to 


T ~ 


1 


87 ra(t)M ’ 


(48) 


which shows that the a(t)M plays the role of mass, and is compatible with the energy accretion in the conformal BHs 
fl9| - {2lj . [ 51 ] . Again, we see that the temperature analysis can support our expectation from M as the physical mass 
in our solutions. For the area of quasi horizon, we have 


A = J \/addd(j) = 4na(tfrl = I67 tM 2 . 


(49) 


In the mcVittie spacetime, this integral leads to A = 16ttAI 2 . In order to vindicate our approximation, we consider 
S = A for the entropy of quasi BH. In continue and from Eq. (S3, we get 


TdS ~ dM = dE. 


(50) 


Whereas, we reach to TdS ~ dM ^ dE for the mcVittie spacetime. In the coordinates (i, rj , 9 , <j>), we should remind 
that, unlike the mcVittie spacetime, E = Mqw = Mms — M is valid for quasi BH and the work term can be 
neglected as the result of slow expansion approximation (dW ~ 0) [5lj] . Finally and unlike the mcVittie’s horizon, 
we see that TdS — dE is valid on the quasi event horizon. This result points us to this fact that the first law of 
the BH thermodynamics on quasi event horizon will be satisfied if we use either the Gong-Wang or the Misner-Sharp 
definitions for the energy of quasi BH. TdS ~ dE is valid for the conformal Schwarzschild BH, too [51]. For the flat 


background, we see that the surface area at redshift singularity in our spacetime is equal to the mcVittie metric which 
is equal to the Schwarzschild metric. In continue and by bearing the slow expansion approximation in mind, we saw 
that the temperature on quasi horizon is like the Schwarzschild spacetime [191 ] . In addition, we saw that the quality 
of the validity of the first law of the BH thermodynamics on quasi event horizon is like the conformal Schwarzschild 
BH’s and differs from the mcVittie’s solution. 








In another approach and for the mcVittie spacetime, if we use the Hawking-Hayward definition of mass as the total 
confined energy to the hypersurface f = , we reach to 


TdS ~ dM = dE, (51) 

where we have considered the slow expansion approximation. In addition, Eq. m will be not valid, if one uses 
the Komar mass (1331) . Finally, we saw that the first law of thermodynamics will be approximately valid in the 
mcVittie’s solution, if one uses the Hawking-Hayward definition of energy. Also, none of the Komar, Misner-Sharp 
and Gong-Wang masses can not satisfy the first law of thermodynamics on the mcVittie’s horizon. 


IV. OTHER POSSIBILITIES 


According to what we have said, it is obvious that there are two other meaningful sentences in expansion (1241) . The 
first term is due to n = — 2 and points to the charge, where the second term comes from n = 2 and it is related to 
the cosmological constant. Therefore, the more general form of f(t,r ) can be written as: 


= 1 - 


2 M 


+ 


Q 2 


a(t)r (a(f)r) 


~ 3 A(a(f)r) 2 , 


(52) 


where we have considered the slow expansion approximation and used these definitions 6_2 = Q 2 and 62 = — ^A. 
Imaginary charge (fr _2 < 0) and the anti De-Sitter (A < 0) solutions are allowed by this scheme, but these possibilities 
are removed by the other parts of physics. Consider Eq. (EU) when A = 0, there are two horizons located at 

r+ = ~and r_ = M 0L' These radiuses are same as the Gao’s flat case [39J. In the low expansion 

regime (a(t) ~ c), these radiuses point to the event and the Coushy horizons, as the Riessner-Nordstorm metric [§j. 
Hence, we refer to them as quasi event and quasi Coushy horizons. The case with Q = 0, M = 0 and A > 0 has 
attractive properties. Because in the low expansion regime (a(t) — c), one can rewrite this case as 


ds 2 w -(1 - ^ 7 f)dt 2 + d \ + rj 2 dn 2 . 

3 (1 - 3 V 2 ) 


(53) 


This is nothing but the De-Sitter spacetime with cosmological constant A, which points to the current acceleration 
era. 


Horizons and temperature 


Different f(t,r ) yield apparent horizons with different locations, and one can use Eqs. © and © in order to find 
the location and the temperature of apparent horizon. For every f(t,r), using the slow expansion regime, we get: 

ds 2 s=s - f(rf)dt 2 + + rf dCt 2 . (54) 

fvn) 

Now, the location of horizons and their surface gravity can be evaluated by using Eq. m■ Their temperature is 
approximately equal to Eq. m , or briefly: 


Ti ~ 


f'(v) I 

4-7T 


where (') is derivative with respect to radii 77 and 77 ^ is the radii of i th horizon. 


(55) 


V. CONCLUSIONS 

We considered the conformal form of the special group of the non-static spherically symmetric metrics, where it 
was assumed that the time dependence of the conformal factor is like as the FRW’s. We saw that the conformal BHs 
can be classified as a special subgroup of these metrics. In order to derive the new solutions of the Einstein equations, 
we have imposed perfect fluid concept and used slow expansion approximation which helps us to clarify the physical 
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meaning of the parameters of metric. Since the Einstein tensor is diagonal, there is no energy accretion and thus 
the backreaction effect is zero. This imply that the energy (mass) should be constant in our solutions. These new 
solutions have similarities with earlier metrics that have been presented by others f37j-j39}. A related metric which is 
similar to the special class of our solutions was introduced by mcVittie [33, [3I|. These similarities are explicit in the 
flat case (temperature and entropy at the redshift singularity), but the differences will be more clear in the non-flat 
case (k ^ 0), and we pointed to the one of them, when we debate the redshift. In addition and in the flat case, we 
tried to clear the some of differences between our solution and the mcVittie’s. We did it by pointing to the behavior of 
the redshift singularity in the various coordinates, the mass notion, and thermodynamics. Meanwhile, when our slow 
expansion approximation is broken then there is no horizon for our solutions. Indeed, these objects can be classified 
as naked singularities which can be considered as alternatives for BHs [1(1 [47( . 

For the our solutions and similar with earlier works [37l - (39[ . the co-moving radiuses of the redshift singularities 
are decreased by the expansion of universe. Also, unlike the previous works 1371 - 1391 . the redshift singularities in our 
solutions are independent of the background curvature. By considering the slow expansion approximation, we were 
able to find out BH’s like behavior of these singularities. We pointed to these objects and their surfaces as quasi BHs 
and the quasi horizons, respectively. In continue, we introduced the apparent horizon for our spacetime which should 
be evaluated by considering the FRW background. 

In order to compare the mcVittie’s solution with our mcVittie’s like solution, we have used the three existing 
definitions of mass including the Komar mass, the Misner-Sharp mass ( Mms) and the Gong-Wang mass ( Mqw)■ We 
saw that the notion of the Komar mass of quasi BH differs from the mcVittie’s solution. Also, in our spacetime, we 
showed that the Mms and Mqw masses yield the same result on the apparent horizon and cover the FRW’s result 
in the limiting situations. In addition, using the slow expansion approximation, we evaluated Mms and Mqw on the 
quasi event horizon of our mcVittie’s like solution, which leads to the same result as the Komar mass. In addition, 
we should express that, the same as the mcVittie spacetime, the energy conditions are not satisfied near the quasi 
horizon. 

In addition, we have proved that, unlike the mcVittie’s solution, the first law of thermodynamics may be satisfied 
on the quasi event horizon of our mcVittie’s like solution, if we use the Komar mass or either Mms or Mqw as the 
confined mass and consider the slow expansion approximation. This result is consistent with previous studies about 
the conformal BHs [5l|, which shows that the thermodynamics of our solutions is similar to the conformal BHs’. In 
order to clarify the mass notion, we think that the full analysis of the Hawking-Hayward mass for our solution is 
needed, which is out of the scope of this letter and should be considered as another work, but our resolution makes 
this feeling that the predictions by either the slow expansion approximation or using the suitable coordinates for 
describing the metric for mass, may be in line with the Hawking-Hayward definition of energy, and have reasonable 
accordance with the Komar, Mms and Mqw masses of our solutions. Indeed, this final remark can be supported 
by the thermodynamics considerations and the no energy accretion condition ( Gt r = 0). Moreover, we think that, in 
dynamic spacetimes, the thermodynamic considerations along as the slow time varying approximation can help us to 
get the reasonable assumptions for energy and thus mass. Finally, we saw that the first law of thermodynamics will 
be approximately valid in the mcVittie’s and our solutions if we use the Hawking-Hayward definition of the mass in 
the mcVittie spacetime and the Komar mass as the physical mass in our solution, respectively. In continue, the more 
general solutions such as the charged quasi BHs and the some of their properties have been addressed. 

Results obtained in this paper may help achieving a better understanding of black holes in a dynamical background. 
From a phenomenological point of view, this issue is important since after all, any local astrophysical object lives in 
an expanding cosmological background. Finally, we tried to explore the concepts of mass, entropy and temperature 
in a dynamic spacetime. 
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